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ABSTRACT 

An analysis is made to investigate the mass transfer over a stretching surface embedded in a porous 
medium in the presence of first order chemical reaction. Using similarity transformation, the 
governing partial differential equations are transformed into a set of ordinary differential equations 
which solved by shooting method. Comparisons with previously published work on special cases of 
the problem are performed and the results are found to he in excellent agreement. It is observed that 
the local mass transfer and concentration profile are very sensitive to change in the values of 
traction rate parameter, permeability parameter and Schmidt number. 
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I. INTRODUCTION 

The boundary layer flows due to stretching surface in porous medium are relevant to many 
engineering problems such as paper production, preparing plastic and metal sheets etc. The dynamics of the 
boundary layer flow over a stretching surface originated from the pioneering work of Sakiadis [1, 2] who 
initiated the study of boundary layer flow over a continuous solid surface moving with constant speed. Crane 
[3] extended it to analyze the steady two dimensional boundary layer flow caused by the stretching of elastic 
flat surface which moves in its plane with velocity varying linearly with distance from a fixed point. Many 
authors presented some mathematical results, and good amount of references can be found in the papers by Ali 
[4] and [5], Elbashbeshy [6], Ishak et al. [7] and Elbashbeshy and Bazid [8] .The studies carried out in these 
papers in the case steady state flow. The unsteady state problem over a stretching surface, which is stretched 
with a velocity that depends on time is considered by Anderson et al. [9], Elbashbeshy and Bazid [10] and 
Ishak et al. [11]. The effects of chemically reactive solute distribution on fluid flow due to a stretching surface 
also bear equal importance in engineering researches. The chemical reaction effects were studied by many 
researchers on several physical aspects. The diffusion of a chemically reactive species in a laminar boundary 
layer flow over a flat plate was demonstrated by Chambre and Young [12]. The effect of transfer of 
chemically reactive species in the laminar flow over a stretching sheet explained by Andersson et al. [13]. 
Takhar et al. [14] analyzed the flow and mass transfer on a stretching sheet with a magnetic field and 
chemically reactive species with n-th order reaction. 

The mass transfer in boundary layer flow due to a stretching surface in porous medium also has 
important applications in many industrial problems. The effect of mass transfer in laminar flow over a 
stretching surface was investigated by Radwan and Elbashbeshy [15]. Radwan and Elbashbeshy [15] analyzed 
the flow and mass transfer on a stretching surface with a magnetic field. Akyildiz et al. [16] reported a solution 
for diffusion of chemically reactive species in a flow of a non- Newtonian fluid over a stretching sheet 
immersed in a porous medium. El- Aziz [17] explained unsteady flow due to a stretching sheet with mass and 
heat transfer. Recently, Krishnendu [18] studied the boundary layer flow with first order chemical reaction 
over a porous flat plate. Krishnendu [19] studied the mass transfer on a continuous flat plate moving in a 
parallel or reversely to a free stream in the presence of a chemical reaction. Ferdows et al. [20] investigated the 
effects of order of chemical reaction on mass transfer over a linearly stretching surface. Based on the above- 
mentioned investigations and applications, this paper is concerned with two-dimensional steady, 
incompressible, laminar boundary layer flow of a fluid over a linearly stretching surface. 

In this paper we investigate numerically the effects of chemical reaction on the steady laminar two- 
dimensional boundary layer flow and mass transfer over a stretching surface embedded in porous medium. 
The method of solutions based on the well-known similarity analysis together with shooting method. 
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II. Formulation of the problem 

Consider two dimensional steady, laminar boundary layer flow of a fluid over a linearly stretching surface 
(i.e. stretched with a velocity proportional to x) embedded in porous medium with velocity u and 
concentration C moving axially through a stationary fluid. The fluid is viscous incompressible with constant 
physical properties. We assume that the fluid far away from the surface is at rest and at concentration Coo. . 
The x - axis runs along the continuous surface in the direction of motion and y - axis is perpendicular to it. 
The continuity, momentum and reactive concentration equations for governing the flow and concentration 
distribution in the boundary layer region along the stretching surface may be written as 

du 8v 

+— =0 (1) 

dx dy 



The boundary conditions are given by 

y = : u = Ax, v = 0, C = C a + Ax " 

(4) 

y->oo:u = 0,C=C B 

where u and v are the velocity components in the x and y directions, respectively, k is the permeability of 
porous medium. C is the concentration of the fluid, v is the kinematic viscosity, R is a constant of first-order 
chemical reaction rate, n is a power-law exponent, which signifies the change of amount of solute in the .re- 
direction, D is the effective diffusion coefficients and A and X are constants. The equation of 
continuity is satisfied if we choose a dimensionless stream function y/ ( x , y ) such that 

dy ' dx' 

In order to solve Equations (l)-(4), we introduce the following similarity transformation and dimensionless 
concentration. 

V(x,y) n C-C m 



4 



where rj is the similarity variable,^ (x, >•) is the dimensionless stream function, 0(tj) is the dimensionless 
concentration. Consequently, equations (2)-(3) and the boundary conditions (4) can be written in the following 
form 

/"' + //"- f' 2 ~ Kf '= (5) 



+ Sc ( fO - nf ) - ScR c = 
tj = : / = 0,/' = 1,0 = 1 



rj -> oo : / ' = 0,0 = 
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Where prime denotes differentiation with respect to tj . K = is a permeability parameter, R = — is 

nk ' n 

reaction rate parameter of the solute and Sc = — is Schmidt parameter. 

D 

The mass transfer coefficient in terms of the Sherwood number Sh is given by 



8C 




where R e = "' is the local Reynolds number. 

III. Numerical method 

The above Eqs. (5) and (6) along with boundary conditions (7) are solved by converting them to an initial 
value problem. We set 



p ' = Kz + z 2 - fp (8) 
q = Sc[R0 + nz0- fq] (9) 
with the boundary conditions 

/(0) = 0, /'(0) = 1, 0(0) = 1 (10) 

In order to integrate (8) and (9) as an initial value problem we require a value for p(0) i.e. /"(0) 
and (9 ' (0) but no such values are given in the boundary. The most important factor of shooting method is to 
choose the appropriate finite values of rj x .In order to determine ?; rj for the boundary value problem stated by 
Eqs. (8)- (9), we start with some initial guess value for some particular set of physical parameters to obtain 
/ " (0) and 6* ' (0) . The solution procedure is repeated with another large value of tj ^ until two successive 
values of / " (0) and 6* ' (0) differ only by the specified significant digit. The last value of tj x is finally chosen 
to be the most appropriate value of the limit rj _ r for that particular set of parameters. The value of rj a may 
change for another set of physical parameters. Once the finite value of rj x is determined then the integration is 
carried out. We compare the calculated values for / ' and at rj = 10 (say) with the given boundary 
conditions / ' (10 ) = and (10 ) = and adjust the estimated values, / " (0) and# ' (0) , to give a better 
approximation for the solution. We take the series of values for / " (0) and# ' (0) , and apply the fourth order 
classical Runge-Kutta method with step-size h= 0.01. The above procedure is repeated until we get the results 
up to the desired degree of accuracy, 10 5 . 
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IV. RESULTS AND DISCUSSION 

The set of non-linear ordinary differential equations (5) and (6) satisfying the boundary conditions (7) 
have been solved numerically using the Mathematica method for several values of the involved parameters, 
namely permeability parameter, power-law exponent, reaction rate parameter of the solute and Schmidt 
parameter. In case n=0 and K = our results are very similar to those considered by Uddin et al. [21], Takhar 
et al. [14] and Andersson et al. [13]. To validate the numerical method used in this study, the case n=0 and 
K = , was considered and the results for concentration gradient are compared with the numerical solution 
which reported in Uddin et al. [21], Takhar et al. [14] and Andersson et al. [13]. The quantitative comparison is 
shown in Table 1, and found to be in excellent agreement. 

Table 1 Comparison of the values of- 6 ' (0) with that of Takhar et al. [14], Andersson et al. [13] and Uddin et 
al. [21] forn=0andK = . 



Present Study 
0.149083 



Andersson et al. [13] 
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Fig. 1: Velocity profiles/ (n) for various values of K 
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Sc = 1 , Rc = 1 , n = 1 



K = 10 , 4 ,1 ,0 



Fig. 2: Concentration profiles (7; ) for various values of K 

The numerical solution of velocity has presented for various values of the permeability parameter K . The 
reactant solute equation is solved numerically and the results are shown graphically. The velocity profiles for 
various values of the permeability parameter K have been plotted in Fig. 1. From the figure it is noted that with 
increase of K , the velocity for any fixed value of n decreases. 

1.2 





K = 1 , Rc = 1 , n = 1 




Sc = 1 , 2 , 5 , 1 1 
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Fig. 3. Concentration profiles (tj ) for various values of Sc 
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Fig4. Concentration profiles (;; ) for various values of R c 

Thus it is clear that the permeability parameter K opposes motion. Consequently, the momentum 
boundary layer thickness reduces with the increase K in and this fact is also seen from wall shear stress 
behaviour. Exhibits concentration profiles for various values of K . The value of contaminate solute at 
particular value of r| increases with the increase of the permeability parameter K and also the concentration 
boundary layer thickness increases. This implies that the permeability parameter acts to enhance the distribution 
of the reaction solute on the stretching surface. This result may be useful, in the situation where the 
enhancement of solute transfer from the surface is the prime important. 



1.2 





K =1 , Sc = 1 , Rc = 1 






n = , 2 , 6 , 10 
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Fig. 5. Concentration profiles (rj ■ various values of w(>0). with K=Sc= R c =1 
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Fig. 6. Concentration profiles 9 (77 ) for various values of «(<0). with K=Sc= R c =1 

Now, we concentrate on variation the solute curves for different values of Schmidt number Sc. The 
curves are drawn in the Fig. 3. The Schmidt number has major effects on the distribution of solute. The 
concentration boundary layer thickness as well as the concentration at a fixed point decreases quickly with 
increasing values of Sc. This is due to the fact that the rate of solute transfer from the surface increases when the 
Schmidt number increases. Figure 4 is the graphical representation of concentration profiles for various values of 
reaction rate parameter R c .It has been found that the reaction rate parameter affect the solute profiles in similar 
way as that of the Schmidt number i.e., the increase of R c reduces both the solute boundary layer thickness and 
value of the solute at fixed r\. So, in case of the distribution of reactive solute, the reaction rate parameter is a 
decelerating agent. 

Finally, Figs. 5 and 6 exhibit the concentration profiles in the boundary layer flow region for different 
values of power-law exponent n. It is noticed from Fig. 5 that for the increasing values of n with n>0, the curve 
representing the distribution of solute for specific value of r| decreases. While, in Fig. 6 the concentration profile 
increases with increase in the magnitude of n with n<0 and for large negative values of n, the overshoot of solute 
is observed near the surface. Thus, the effect of increase of n when the surface concentration is 
C w = C m + Ax " is completely opposite to the effect of increase n when the surface concentration 
is C w = C r + A I x " where A is constant positive. Note that, the wall concentration is constant when n=0. 

Table 2 

Values of skin-friction coefficient / " ( ) and - 6 1 ( ) for various values of permeability parameter K . 



K 


/"(0) 


6>'(0) - 





-1.00000 


1.45942 


1 


-1.414214 


1.414214 


4 


-2.236068 


1..3433 


10 


-3.316625 


1..2781 
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Table 3: Values of - (0 ) for various values of Schmidt number Sc 



Sc 


0'(O)- 


1 


1.414214 


2 


2.063933 


5 


3.. 366092 


11 


5.085476 



Table 4: Values of - (0 ) for various values of reaction rate parameter R c 



R c 


0'(O)- 





0.892503 


1 


1.414214 


4 


2.2692285 


12 


3.641534 


Table 5: Values of - 9 ' (0) for various values of n 


n 


0'(O)- 





1.154956 


2 


1.649948 


6 


2.434269 


10 


3.061963 



Tables2-5summarize the calculated values of local skin friction coefficient / " (0) and rate of mass transfer - 
f?'(0) for 1 < 5c < 11 , < n < 10 and < R c < 12 . These tables indicate that - /" (0) increasing as 
permeability parameter K increases. Since K does not occur explicitly in the diffusion equation , it effect on 
mass transfer -6 1 (0) - is small. A lso, these tables indicate that - 6 ' (0) increasing as n, R c and Sc increases. 



V. CONCLUSIONS 

In this investigation, an analysis is made to find the behaviour of the distribution of reactive solute 
undergo a first order reaction in steady boundary layer flow incompressible fluid over a stretching surface 
embedded in porous medium with variable surface concentration. Using the similarity transformation a set of 
ordinary differential equations is obtained from the governing equations. The momentum equation and the 
equation of reactive solute are solved numerically. The results show that the permeability parameter tends to 
reduce the rate of flow from the wall and is broadening the solute layer. The Schmidt number and the reaction 
rate parameter reduce the solute boundary layer thickness. Most, importantly, the effects of initial variable 
solute distribution over a stretching surface is interesting i.e. for the increase in magnitude of n, the 
concentration decreases when «>0 whereas increases when n<0. 
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